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HETEROCLINIC CYCLES ARISING IN GENERIC UNFOLDINGS OF NILPOTENT 

SINGULARITIES 

PABLO G. BARRIENTOS, SANTIAGO IBANEZ, AND J. ANGEL RODRIGUEZ 


Abstract. In this paper we study the existence of heteroclinic cycles in generic unfoldings 
of nilpotent singularities. Namely we prove that any nilpotent singularity of codimension 
four in unfolds generically a bifurcation hypersurface of bifocal homoclinic orbits, that 
is, homoclinic orbits to equilibrium points with two pairs of complex eigenvalues. We also 
prove that any nilpotent singularity of codimension three in R^ imfolds generically a bifurca¬ 
tion curve of heteroclinic cycles between two saddle-focus equilibrium points with different 
stability indexes. Under generic assumptions these cycles imply the existence of homoclinic 
bifurcations. Homoclinic orbits to equilibrium points with complex eigenvalues are the sim¬ 
plest configurations which can explain the existence of complex dynamics as, for instance, 
strange attractors. The proof of the arising of these dynamics from a singularity is a very 
useful tool, particularly for applications. 

1. Introduction 

The relationship between dynamic complexity and the presence of homoclinic orbits was 
discovered by Poincare more than a century ago. In his famous essay on the stability of 
the solar system [45], Poincare showed that the invariant manifolds of a hyperbolic fixed 
point of a diffeomorphism could cut each other at points, called homoclinics, which yield 
the existence of more and more points of this type and consequently, a very complicated 
configuration of the manifolds. Many years later, Birkhoff [6] showed that, in general, near 
a homoclmic point there exists an extremely intrmcated set of periodic orbits, mostly with 
a very long period. By the mid 60's, Smale [54] placed his geometrical device, the Smale 
horseshoe, in a neighborhood of a transversal homoclinic point. The horseshoes explained 
the Birkhoff's result and arranged the complicated dynamics that occur near a homoclinic 
orbit by means of a conjugation to the Bernouilli's shift. In [40] authors proved the appearance 
of strange attractors during the process of creation or destruction of the Smale horseshoes 
which appear through a bifurcation of a tangential homoclmic point. These attractors are 
like those shown in [2] for the Henon family, that is, they are nonhyperbolic and persistent 
in the sense of measure. 

In the framework of vector fields, Shil'nikov [52] proved that in every neighborhood of a 
homoclinic orbit to a hyperbolic equilibrium point of an analytical vector field on R^, with 
eigenvalues A and -g ± coi such that 0 < p < A, that is, the so-called Shil'nikov homoclinic orbit, 
there exists a countable set of periodic orbits. This result is similar to that found by Birkhoff 
for diffeomorphisms and thus, it should be understood in a maimer similar to that devised 
by Smale. Indeed, Tresser [55] showed that in every neighborhood of such a homoclinic 
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orbit, an infinity of linked horseshoes can be defined in such a way that the dynamics is 
conjugated to a subshift of finite type on an infinite number of symbols. Once again, these 
horseshoes appear and disappear by means of generic homoclinic bifurcations leading to 
persistent non hyperbolic strange attractors like those in [40]. 

As follows from [41], nonhyperbolic dynamics is dense in the space X of vector fields with 
a Shil'nikov homoclinic orbit. In particular, for each £ > 0, the subset of vector fields with a 
homoclinic tangency to a hyperbolic periodic orbit in an £-neighbourhood of the homoclinic 
orbit is dense in X. These tangencies give rise to suspended Henon-like strange attractors. 
In [46, 47] it was proved that infinitely many of these strange attractors can coexist in non 
generic families of vector fields with a Shil'nikov homoclinic orbit, for parameter values 
in a set of positive Lebesgue measure. Later [25], it was proved that an infinity of such 
attractors can coexist in a more general context. For an extensive study of the phenomena 
accompanying homoclinic bifurcations, see [4, 26, 48]. 

Because of the importance of homoclinic orbits in Dynamics, many papers were devoted 
to prove their existence. A seminal work was due to Melnikov [37], who introduced original 
ideas to prove the existence of transversal homoclinic orbits in non-autonomous perturba¬ 
tions of a planar hamiltonian vector field. These ideas were developed in [10] in order to 
determine both, homoclinic bifurcation curves and the existence of subharmonics in two- 
parameter families of non-autonomous second order differential equations. In [43], Palmer 
developed a theory involving transversal homoclinic points and exponential dichotomies 
that was very useful for the study of homoclinic bifurcations in higher dimensions. 

Since Shil'nikov homoclinic orbits are not transversal, Melnikov's techniques had to be 
modified in order to prove their existence in families of vector fields. In [49], generic families 
of quadratic three dimensional vector fields with Shil'nikov homoclinic orbits were given. 
Putting together ideas from [49, 10, 43], it was proved in [28] that Shil'nikov homoclinic 
orbits appear in generic unfoldings of a nilpotent singularity of codimension four in . Since 
singularities are the simplest elements to be found in phase portraits of vector fields, arguing 
the existence of homoclinic orbits from the presence of singularities is a highly relevant task. 
Nevertheless, in order to get the greatest interest in applications, such singularities should 
be of codimension as low as possible. With this in mind, the result obtained in [28] was 
improved in [29], where it was showed that 

Theorem A. Shil'nikov homoclinic orbits appear in every generic unfolding of the nilpotent singu¬ 
larity of codimension three in R^. 

Proving that Shil'nikov homoclinic orbits can be unfolded generically from a singularity of 
codimension less than three is currently a very interesting open problem. The dimension of 
the center manifold should be at least three. The lowest codimension singularities in R^ with 
a 3-dimensional center manifold are the Hopf-zero singularities which have codimension 
two [22]. The difficulties that appear on studying the existence of Shil'nikov homoclinic 
orbits in generic unfoldings of Hopf-zero singularities are discussed in [17]. 
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Theorem A was essential in [18] to prove the existence of persistent strange attractors in 
the four parametric family of vector fields obtained when two Brusselators are linearly cou¬ 
pled by diffusion. Indeed, this family is a generic unfolding of three-dimensional nilpotent 
singularities of codimension three. Therefore it displays ShiTnikov homoclinic orbits and, 
consequently, persistent strange attractors. Nevertheless, this family may display a richer 
dynamics. Three-dimensional nilpotent singularities appear along two bifurcation curves 
which emerge from a bifurcation point corresponding to a four-dimensional nilpotent sin¬ 
gularity of codimension four, for which the family is also a generic unfolding. Therefore, 
one should wonder whether a different class of homoclinic orbits can take place from this 
four-dimensional nilpotent singularity. In this paper we will prove the following result: 

Theorem B. In every generic unfolding of a four-dimensional nilpotent singularity of codimension 
four there is a bifurcation hypersurface of homoclinic orbits to equilibrium points with two pairs of 
eigenvalues pk ± coki, with k = 1,2, such that pi < 0 < p 2 - 

Homoclinic orbits in Theorem B are usually known as bifocal homoclinic orbits or, shortly, 
bifocus. ShiTnikov [53] was again the first one in studying the dynamics associated to them. 
He proved, as in [52], the existence of a countable set of periodic orbits in the non-resonant 
case -pi ^ p 2 . Subsequent works [13, 21, 33, 23] were devoted to analyze the formation and 
bifurcations of these periodic orbits by studying the Poincare map associated to the flow 
in a neighborhood of the bifocus. Devaney [13] considers the hamiltonian case, hence with 
“Pi = Pi- He proves that for any local transverse section to the homoclinic orbit, and for 
any positive integer N, there is a compact invariant hyperbolic set on which the Poincare 
map is conjugate to the Bernoulli shift on N symbols. In seeking to determine the invariant 
set of this Poincare map in the general case, it is shown in [21] that this set is contained 
in a neighborhood of a spiral sheet (shaped like a scroll). In fact, the invariant set is a 
neighborhood of the intersection of this scroll and its image under the map, which is another 
scroll, in general skewed and offeet from the original. In [33] the authors extend the known 
theory regarding bifocal homoclinic bifurcations and present numerical verification of the 
more interesting theoretical predictions that had been made. Harterich [23] studies bifocal 
homoclinic orbits arising in reversible systems, hence again with-pi = p 2 . He proves that for 
any N >2 there exists infinitely many N-homoclinic orbits in a neighborhood of the primary 
homoclinic orbit. Each of them is accumulated by one or more families of N-periodic orbits. 

As for ShiTnikov homoclinic orbits, it has been proved (see [42]) that homoclinic tangencies 
to hyperbolic periodic orbits are dense in the space of vector fields with a bifocal homoclinic 
orbit. Nevertheless, despite the abundant literature regarding bifocus, as far as we know, no 
result has been established relating the existence of bifocal homoclmic bifurcations with the 
existence of persistent strange attractors. This, in spite of a bifocus seems to be a scenario for 
more complicated dynamics than those inherent to ShiTnikov homoclmic orbits, where the 
existence of such strange attractors has been proved. In fact it seems natural to think that 
the dynamical complexity associated with homoclinic cycles increases with dimension. For 
instance, strange attractors with more than one positive Lyapunov exponent could appear. 
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Therefore, bearing in mind a possible extension of Theorem B to higher dimensions, we will 
begin its proof working with n-dimensional nilpotent singularities of codimension n. 

Let Xbe a C°° vector field on R” with X(0) = 0 and 1-jet linearly conjugated to L^=i ^k+i^ldx]^. 
Introducing appropriate coordinates, X can be written as 


n-l 

L 

k=l 


. 


with/(x) = 0(1 |x||^) where X = (xi,...,x„). ItissaidthatXisanilpotentsingidarity of codimension 
n if the generic condition d^fjdx^ 0 is fulfilled. As we will explain in Section 3, in 
appropriate coordinates and after rescaling, any generic n-parametric unfolding Xa of a 
nilpotent singularity can be written in a neighborhood of the origin as 


n-l 


k=l 


L +<’'*+ L’'‘» +yh o(‘))^. 


k=l 


( 1 . 1 ) 


with £ > 0 and -I-... -I- = 1. The limit family obtained for £ = 0, will play a main role. It 

is time reversible with respect to the involution R(i/i, y 2 , ■■■, yn) = ((“l)”i/i/ (“l)”~^J/ 2 / -yn) 
for parameter values on the set 

T = {(vi,...,v„) € : v„_ 2 ,- = 0 with / = 0,..., L(n - 2)/2j}, 

where [.J denotes the floor function. This manifold 7“ of dimension [n/2 - IJ is called the 
reversibility set of the n-dimensional nilpotent limit family. 


For n = 4 and for values of the parameters vi < 0, V 2 = V 4 = 0, the limit family can be 
transformed in 

a a a ^ 


d d d , 2 ^ ^ 


( 1 . 2 ) 


where 7]3 = 2 “^^^(-vi)“^^^V 3 . Denoting u = xi and P = -rjs the vector field (1.2) is equivalent 
to the fourth-order equation -I- Pu”{t) + u{t) — u^{t) = 0, which has been widely studied 

[1, 12, 3, 8 ] due to its role in some applications as the study of travelling waves of the 
Korteweg- de Vries equation 

Wt = Uxxxx ~ bU-xxx + 2wMx, 


or the description of the displacement of a compressed strut with bending softness resting on 
a nonlinear elastic foundation [9]. In particular, according to [1], when ?]3 = 2 the vector field 
( 1 . 2 ) has a homoclinic orbit to a hyperbolic equilibrium point at which the linear part has a 
pair of double real eigenvalues +1. We will complete the proof of Theorem B by proving that 
the homoclinic orbit persists for parameter values on a hypersurface TTom which intersects 
£ > 0 and is obtained by studying the appropriate bifurcation equation. Moreover we will 
show that 'Horn contains regions corresponding to bifocal homoclinic orbits. An essential 
fact used in [ 1 ] to prove the existence of homoclinic orbits in ( 1 . 2 ) is that it is a family of 
hamiltonian vector fields. This permits to apply the general theory developed in [24]. Again 
bearing in mind a possible extension of Theorem B to higher dimensions we will prove that, 
for any even n and for parameter values on the reversibility set, the vector fields in the limit 
family of ( 1 . 1 ) are hamiltonian. 
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Methods used in the proof of Theorem B also allow us to prove the existence of topological 
Bykov cycles, which will be defined below, in the case n = 3, according with the following 
result. 

Theorem C. In every generic unfolding of a three-dimensional nilpotent singularity of codimension 
three there is a bifurcation curve of topological Bykov cycles. 


For n = 3 and when vi < 0 and V 2 < 0, the family (1.1) can be transformed into the family 

where = Ivijv^ and v = V 3 /(-l/v 2 )^^^. When v = 0 this limit family is equivalent to the 
third-order equation x"'{t) -I- x'{t) + x{t)^l2 = c^, which has been studied extensively in the 
literature (see [31, 32, 38] and references therein) because it plays a very relevant role in the 
study of the existence of steady solutions and travelling waves of the Kuramoto-Shivashinsky 
Ut + Uxxx + Uxx + w^/2 = 0. For each value of c, family (1.3) has two saddle-focus equilibria P+. 
In particular, the eigenvalues A and -p ±ioj at P- satisfy that 0 < p < A , that is, the spectral 
assumptions in ShiTnikov's theorem. In [31] it is proved the existence of a heteroclinic 
cormection Fi = W“(P_) n W®(P+) when c = Cjt = 15 ^f22|l^. For this same value of c, it 
was proved in [29] the existence of a topologically transverse intersection F 2 between the 
two-dimensional invariant manifolds W“(P+) and W®(P-). Therefore, T = Ti U T 2 U (P_,P+} 
is a heteroclinic cycle, known in the literature as Bykov cycle. Bifurcations arising from the 
breaking of this codimension two heteroclinic cycle have been widely studied in the literature 
[15,20]. In particular, the birth of homoclinic orbits was considered in [15] for the case which 
concerns us. It follows from these papers that in any generic 2-parametric unfolding of 
a Bykov cycle there exist two homoclinic bifurcation curves, which spiral in to the Bykov 
cycle bifurcation point. Genericity can be read in terms of the transversality between the 
2 -dimensional invariant manifolds, which can be replaced by the condition of topological 
transversality, and the generic splitting of the connection along the 1 -dimensional invariant 
manifolds (details can be seen in [15, 20, 29]). 


In this paper we will prove that family (1.3) exhibits topological Bykov cycles for pa¬ 
rameters along a curve y = ((c(£),v(£), c) : e € [0 , Cq)} for some £0 > 0 with c(0) = c^ and 
v(0) = 0. Fixing a section Z, transverse to Fi we consider the splitting function h{c,v,e) = 
{hi{c,v, e),h 2 {c,v,e)) defined as the distance between W“(P_(t)) n Z and W®(P+(t)) n Z for 
T = (c,v, f) close enough to = (C;t/0/0)- Note that h take values on and existence of 
Bykov cycle is equivalent to h{c,v,e) = (0,0). We will prove that 


dhi 


in) 


dhi 


dh2 . , dh2. , 


^ 0 . 


(1.4) 


Hence the existence of y, and therefore Theorem C, follows from the Implicit Function 
Theorem. 
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We should remark that the above generic condition guarantees the existence of Shil'nikov 
homoclinic orbits, in the sense of Theorem A. Indeed, for each £ > 0 small enough and fixed, 
(1.4) implies that the splitting function is a local diffeomorphism or, in other words, that the 
splitting of the cormection along the 1-dimensional invariant manifolds is generic. Therefore 
there exists a Shil'nikov bifurcation surface shaped as a scroll around y. In the proof of 
Theorem A given in [29], for the sake of brevity, we did not include the computation of the 
above generic condition, although we appealed to it. Since the publication of that paper 
such computation has been frequently demanded to us. 

The paper is organized as follows. In §2 we include a brief summary of results about 
dichotomies in order to get a precise formulation of the bifurcation equations which are 
required in the subsequent sections to prove the existence of heteroclinic and homoclinic 
orbits. §3 is devoted to introduce nilpotent singularities on R” and the limit families obtained 
after rescaling properly a generic unfolding. There we also state that when n is an even 
number and the parameter values belong to the reversibility set, the limit family consists of 
hamiltonian vector fields. Theorem B and Theorem C are proved in §4 and §5, respectively. 

2. Dichotomies and bifurcation equations. 

Let x' = f{x) be a nonlinear equation, where x e R” and / is a regular enough vector 
field, and assume that it has a heteroclinic orbit y = {p{t) : f € R] cormecting two hyperbolic 
equilibrium points p+ and p- (if p+ = p-, y is said homoclinic). Consider a family 

x'= f{x) +g{A,x), (2.1) 

with A e R^ and g regular enough, such that ^(0, x) = 0. For any A small enough, family (2.1) 
has hyperbolic equilibrium points p+(A) and p-(A), continuation of p+ and p-, respectively, 
and the stability index is preserved. In order to study the persistence of the heteroclinic orbit 
for A small enough we introduce the change of variables x{t) = z{t) -l- p{t) in (2.1) to obtain 

z'(f) = D/(p(f))z(f) -I- &(A,f,z(f)), (2.2) 

where 

b{A, t, z(f)) = fipit) + z(f)) - f{p{t)) - Df(p(t))z(t) + g{A, pit) + z(f)). 

Notice that b(0, t, 0) = Dzb{0, f, 0) = 0 for all f e R. 

Persistence of heteroclinic orbits in (2.1) implies the existence of bounded solutions for 
(2.2) which, in turn, implies the existence of bounded solutions for a equation as 

z'(t) = Df(p(t))z(t) + bit), (2.3) 

where b belongs to the space C°(R,R”). In the sequel C^(R, R”) denotes the Banach space 
of bounded continuous R'^-valued functions whose derivatives up to order k exist and 
are bounded and continuous. The existence of bounded solutions of a linear equation 
x' = Ait)x + b{t), as that in (2.3), will be given in terms of exponential dichotomies of 
the homogeneous equation x' = Ait)x and its adjoint w' = —AHyw, where A{ty denotes the 
conjugate transpose of A(f). The classical references for the study of exponential dichotomies 
are [39,11,43,44]. 
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2.1. Exponential dichotomy. Let X(t) be a fundamental matrix of 

x' = A{t)x, X e R", (2.4) 

where A{t) is defined and continuous on an interval / c R. 

Definition 2.1. It is said that the equation (2.4) has an exponential dichotomy on J if there exists a 
projection P : R" —> R”, that is, annby n matrix P with = P, and positive constants K, L, a and 
(I such that for every s,t e J, 

||X(f)PX-i(s)|| < for t > s, 

(2 5) 

||X(f)(J - P)X-i(s)|| < for s>t. 

Let us define Ji^{s) = X(s)PX~^(s) for each s € /. Notice that, according with the above 
definition, I^(s) is the projection corresponding to the fundamental matrix Y(t) = X(t)X~^(s) 
of (2.4) and we can give an alternative definition of exponential dichotomy. 

Definition 2.2. It is said that the equation (2.4) has an exponential dichotomy on J if for all s e J 
there exists a projection lP(s} : R” ^ R" and positive constants K, L, a and (I independents ofs such 
that for all t e J the matrix X~^{t)£P{t)X{t) has constant coefficients and 

||X(f)X-i(s).?^(s)|| < for all t > s, 

||X(f)X-i(s)(J - ^{s))\\ < for all s > t. 

Although the notion of exponential dichotomy is stated for any / c R, the most interesting 
cases are when / is not bounded. We are particularly interested in / = [t, oo) or / = (-oo, t]. 
In such cases the notions of stable and unstable subspaces can be introduced in terms of the 
ranges of the projections of the exponential dichotomies. 

Definition 2.3. Suppose that the matrix A{t) in (2.4) is defined and continuous on J = [t, oo) (resp. 
J = {-oo, t]). For each to e J the stable (resp. unstable) subspace for initial time t = to is defined as 

= {^ e R'^ : ||X(f)X-i(fo).^|| ^ 0 when t ^ c«} 

{resp. El = {^ e R'^ : ||X(f)X-i(fo).^|| ^ 0 when t ^ -oo}). 

Below we give a collection of results which can be helpful to follow the paper. Their proofs 
are available in the literature. 

Proposition 2.4. Suppose that the equation x' - A{t)x has an exponential dichotomy on J. 

i) When J = [t, oo), coincides with the range 'R{SP{to)) of SP{to) for all to e J. Furthermore 

n^to)) = {^ € R" : sup ||X(f)X-i(fo)-^ll < 

t>to 

and for all fo, L c / it follows that E^ = X{ti)X~^{to)El. 

ii) When ] = {-oo, z], El coincides with the kernel N{l3^{to)) of SP{to) for all to e /. Furthermore 

N{lP{to)) = {Ee R” : sup ||X(f)X-i(fo)^|| < «.}, 

t<to 

and for all to,h £ J it follows that El = X{ti)X~^{to)El. 


BARRIENTOS, IBANEZ, AND RODRIGUEZ 


From the above proposition it follows that the linear flow sends and to and , 
respectively. Accordingly, once and £“ are fixed, the stable and unstable subspaces are 
determined for all t. Therefore, the projections are also determined for each t e J once they 
are defined for f = fo- The same observation follows taking into account the uniqueness of 
solutions for the equation 

^'{s) = X'(s)PX-i(s) + X(s)P(X-Hs))' = A{s)^is) - ^{s)A{s). 

Lemma 2.5. If the linear homogeneous equation x' - A(t)x, with t e (-00,00), has exponential 
dichotomy [t, 00) (resp. (-00, T])/or some t e R then it has exponential dichotomy on [fo, 00) (resp. 
(-00, to]) for all to e R. 


The next result [44, Lemma 7.4] states that exponential dichotomy is a robust property 
with respect to small enough perturbations of A{t). 


Proposition 2.6. Suppose that x' = A{t)x has an exponential dichotomy on } = [a,b] (with -00 < 
a <b < 00 ) with projection matrix function I^{t), with constants Ki, K 2 and exponents ai, a 2 - Let fli 
and II 2 be such that 0 < fli < aiandO < ^2 < « 2 - Then there exists 60 = 5o(Xi,X2/«i/«2,iSi,|32) > 0 
such that ifB{t) is a continuous matrix function with 

||B(f)|| < 6 t < 60 for all t e }, 


the perturbed system 


x' = [A(f) + B{t)]x 


has an exponential dichotomy on J with constants £ 1 , £2 exponents fli, ^2 o-nd projection matrix ^(f) 
satisfying that 

\m) - nm < N6t, 


where £ 1 , £2 and N are constants which only depend on K^, K 2 , ai and a 2 - 


From the above result and Lemma 2.5 it follows the existence of an exponential dichotomy 
for the homogeneous part z' = Df(p{t))z of the equation (2.3). Since 

lim p{t) = p+ and lim p{t) = p_ 

and according to Proposition 2.6, the equation x' = Df{p{t))x has the same exponential 
dichotomy than x' = Df(p+)x (resp. x' = Df{p-)x) on [fo/O®) (resp. (- 00 , tg]). That is, if the 
stable (resp. unstable) subspace of x' = Df{p+)x (resp. x' = Df{p-)x) has dimension k then 
x' = Df{p{t))x has an exponential dichotomy on [to, 00 ) (resp. (- 00 , tg]) with stable subspace 
(resp. unstable subspace £“^) with dimension k. In fact we have the following result: 

Proposition 2.7. Let p{t) he a solution of the equation x' - f{x) parametrizing an orbit on the stable 
(resp. unstable) manifold of an equilibrium point p. Hence the variational equation x' = Df{p{t))x 
has exponential dichotomy on [to, 00 ) (resp. (- 00 , tg]). Moreover, 

'Rmto)) = Tp^t^)Wfp) (resp. Af(^(fo)) = Tp(f„)W''(p)). 
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Now we can apply to (2.3) the result below, which relates the existence of bounded 
solutions for a linear equation and for its adjoint. 


Theorem 2.8. [43, Lemma 4.2] Let A{t) be a bounded and continuous matrix defined on (-oo, oo). 
The linear equation x' = A{t)x has exponential dichotomy on [fo/°o) ^^nd on (-oo, tg] if and only if 
the linear operator 

L : x{t) € C^(]R,R") ^ x'{t)-A{t)x{t) e C°(R,R") 
is Fredholm. The index ofL is dimE®^ + dim£“ - n. Moreover, b e 'R{L) if and only if 



< zv{t),b{t) > dt = 0 


for all bounded solutions zv{t) of the adjoint equation w' = -A{tyzv. 


To explore the existence of bounded solutions of the adjoint equation one has to study its 
properties of exponential dichotomy. 


2.2. Exponential dichotomy for the adjoint equation. Let X{t) be a fundamental matrix of 
the equation x' = A{t)x. It is well known that the conjugate transpose of its inverse X~^{ty is 
a fundamental matrix of the adjoint equation w' = —A{tyw. From this relationship between 
the fundamental matrices of both equations we can conclude the following result about the 
connection between their respective dichotomies. 

Proposition 2.9. If the linear equation x' - A{t)x has exponential dichotomy on J with projection 
matrix SA{t) then the adjoint equation w' = -Ajtfw has exponential dichotomy on J with projection 
matrix I - £P{ty. Moreover, for each to e J 

R" = Kmto)) T m - ntoT) = t Nmkn 
R" = m - T n^ito)*) = Af(^(fo)) ^ 'Rmto)*). 

As done in Definition 2.3 we can define now the stable and unstable subspaces for adjoint 
equations. 

Definition 2.10. Suppose that J = [t, oo) (resp. J = (-oo, t]) is contained in the interval of definition 
ofx' = A{t)x. For each to & ] the stable (resp. unstable) subspace for initial time t = to of the adjoint 
equation x' - -Ajtfx is defined as 

Ejl = {w e R” : ||X“^(f)*X(to)*a;|| ^ 0 when t ooj 
(resp. E'j* = {re e R” : ||X“^(f)*X(fo)*^p|| — > 0 when t —> -ooj). 

The following result about the relationship between the invariant subspaces of the equa¬ 
tion x' = A{t)x and its adjoint follows as a straight consequence of Proposition 2.4 and 
Proposition 2.9. 
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Proposition 2.11. Suppose that the equation x' = A{t)x with a e R” and t e } has exponential 
dichotomy in ]. 

(1) IfJ = [to/Oo) then 

= {A e R'^ : sup ||X(i)X-i(io)^ll < 

t>to 

Ell = = {a; e R” : sup ||X-i(i)*X(io)*a;|| < oo}, 

t>to 

and W = El ± £f. 

to to 

(2) If} = (- 00 , to] then 

El = Umto)) = {A e R^ sup ||X(t)X-i(to)x|| < oo}, 

t<to 

Ef* = n^itoT) = {a; € R" : sup ||X-i(t)*X(to)*a;|| < oo}, 

t<to 

and R" = £“ A . 

to to 

In short, if the linear equation x' = A{t)x has exponential dichotomy in / = [to, oo) (resp. 
(-00, tg]) then the forward (resp. backward) bounded solutions of this equation and its 
adjoint are those which tend to zero exponentially when t —> oo (resp. t —> -oo). On 
the other hand, from the decompositions of R" given in Proposition 2.11 it follows that, if 
a' = A{t)x has m linearly independent forward (resp. backward) bounded solutions, then 
the adjoint equation w' = —A{t)*w has n-m linearly independent forward (resp. backward) 
bounded solutions. 

Proposition 2.12. If the linear equation x' = A{t)x has exponential dichotomy in [fo/^o) (tnd 
in {-oo, fg] then the number of linearly independent bounded solutions of the adjoint equation 
w' = -A(t)*w is 

dimEf n Ef = n — dim£? - dimEj' + dim£? n E't. 

iQ iQ to to ro to 

Now we apply the above result to determine the number of bounded solutions of the 
adjoint equation z' = -D/(p(f))*z. As we have already noticed, the number of linearly 
independent forward (resp. backward) bounded solutions of the variational equation x' = 
Df{p{t))x is given by the dimension of the stable (resp. unstable) subspace of the equation 
a' = Df{p+)x (resp. x' = Df{p-)x). That is, such number coincides with the dimension 
of W^{p+) (resp. W“{p-)). Therefore, taking into account that ^to = TpitomPr) and El = 
Tp(fg)W"(p-), we can conclude, from Proposition 2.12, the following result. 

Proposition 2.13. Ifp{t) is a (homo)heteroclinic solution connecting two equilibrium points p+ and 
P- then the number of linearly independent bounded solutions of the adjoint variational equation 
w' = -Df{p{t))*w is the codimension ofTp(^tg)Wfp+) + that is, 

n - dimW"(p+) - dim W"(p_) + dimTp(t(,)W"(p+) n Tp(t(,)W''(p_). 
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Definition 2.14. A (homo)heteroclinic orbit y is said non degenerate if 

dim TpW^(p+) n TpW“(p_) = 1, 
with p e y. Otherwise y is said degenerate. 

Remark 2.15. If the (homo)heteroclinic orbit is non degenerate, the number of linearly inde¬ 
pendent bounded solutions is obtained directly from the stability indexes of p+ and p_. More¬ 
over, although dimTp(fg)W'®(p+) = dim W®(p+) and dimTp(tg)W"(p-) = dim W''(p-), in general 
dimTp(fg)W®(p+) n Tp(f(,)W"(p-) does not coincide with dim W®(p+) n W“(p-). 

In the sequel the (homo)heteroclinic orbit y = {p{t) : t e R} will be non degenerate. 


2.3. Bifurcation equation. As already mentioned, the existence of (homo)heteroclinic orbits 
for (2.1) implies the existence of bounded solutions of (2.2) and, consequently, the existence of 
bounded solutions of (2.3) when b{t) e C°(R, R”). According to Proposition 2.8, if the adjoint 
variational equation w' = -Df{p{t))*w has d linearly independent bounded solutions Wi, then 
the persistence of the (homo)heteroclinic orbit requires the fulfillment of the d conditions 



{wi{t),b{t)) dt = 0 


for i = 1 ,... ,d. 


The question now is the sufficiency of such conditions. 

When d = 1 the sufficiency could be followed from [10]. In general, for d > 1, the 
techniques to be used follow the first steps of the Lin's method [34, 50]. For ||A|| small 
enough, one has to look for solutions pjj^(-) and p^(-) of (2.1), contained in the stable and 
unstable invariant manifolds of the equilibrium points p+(A) and p-(A), respectively (see 
Figure 1). Initial values p*(fo) will belong to a section transverse to the (homo)heteroclinic 
orbit y. Namely 

= Pito) + imto))}^ = Pik) + (w+ © W- © El) 
where = E®* n E^^ and Wj|j (resp. W“) is the orthogonal complement of 


^ = span(/(p(fo))| in E^ (resp. £“). 


Moreover the condition E°°iA) = p^(fo) “ P^ik) ^ P‘1 will be required. Under these assump¬ 
tions there will exist two unique solutions p*(-) for each A. The jump 


nA) = p-,ik)-ptik) 


measures the displacement between the stable and unstable invariant manifolds on the 
section Ltg in the direction of the subspace E*^ = [E®^ + E“ ]-*-. 

The proof of the result below can be found in [50, Lemma 3.3] and [30, Lemma 2.1.2]. 
Namely, in [30] only the first item is proved and, moreover, the proof is developed for the 
degenerate case although the non degenerate one follows in a similar manner. The second 
item is proved in [50] for the non degenerate case. We include in Appendix C a complete 
and simplified proof of this result. 
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Figure 1. Non-degenerate heteroclinic orbit in where the 1-dimensional manifolds coincide. In 
this case, E® = E'f = Et„ (unidimensional), Ef = E“f = (bidimensional) and = p(fo) + Ef. 
For simplicity, we have assume that the perturbation satisfies g(A,p±) = Ofor all A. 


Lemma 2.16. There exists 6 > 0 such that for all A € with ||A|| < 6 , 

(1) There exists a unique pair of solutions p'^(t) and pf{t) of (2.1) parametrizing orbits on 
W^(p+(A)) and W''(p_(A)), respectively, such that p*(io) ^ ^to 

nA)=pfito)-pt{to)eEl 

Writing the solutions as p* (t) = p(t) +z* (t), then z* (•) are, respectively, forward and backward 
bounded solutions of the equation (2.2). They depend regularly on A and the functions Zp are 
identically zero. 

(2) For e > 0 small enough, there exists a (homo)heteroclinic solution Pa(0 such that 

WpAik) - pito)\\ < £ if and only = 0, 

that is, the components ^“(A) of in a basis {wi: i = 1 ... d} ofE*^ satisfy 

< Wi(s),b(A,s,zf(s)) > ds + I < Wi(s),b(A,s,z^(s)) > ds = 0 

OO %JtQ 

being wfs) = X~^{syX{toywi for i = 1 ,... ,d bounded linearly independent solutions of the 
adjoint variational equation. 

According with the above statement the persistence of (homo)heteroclinic orbits follows 
from the analysis of the bifurcation equation E,°°iA) = 0. The existence of non zero parameter 
values A e R^ such that (A) = 0 follows from the Implicit Function Theorem when Da 
has rank d <k. Thus, the following result follows: 








HETEROCLINIC CYCLES ARISING IN GENERIC UNEOLDINGS OE NILPOTENT SINGULARITIES 


13 


Theorem 2.17. Let = 0, with A e be the bifurcation equation of the differential equa¬ 

tion (2.1). If k > d and rankDA^^CO) = d, then (2.1) has a (homo)heteroclinic orbit for each 
parameter value A on a regular manifold of dimension k-d with tangent subspace at A = 0 given by 
the solutions of the system 

k 

XiJAj-O i = l. d 

;=i 


where 


for i = l,...,d and j 


" f ^ m(s),DAjg(0,pis)) > ds 

kJ — OO 

l,...,k. 


Remark 2.18. Note that, when k < d, A = 0 is the unique value of A e for which there 
exists a (homo)heteroclinic orbit y,\ = {pA{t) : = /(Pa(0) +S'('^/Pa( 0) ^ e R} such that 

supfgj^ |IPa( 0 - p(OII small enough. Ifk>dthe homoclinic connection persists for parameter 
values on a manifold of codimension d where 

d = n - dim W®(p+) - dim W“(p-) + 1. 


In such a case we say that there is (homo)heteroclinic bifurcation of a non degenerate orbit at A = 0 
which is of codimension d. 


3. Nilpotent singularities of codimension n on R” 


3.1. Generic unfoldings. Let X be a C°° vector field in R” with X(0) = 0 and 1-jet at the 
origin linearly conjugated to L^=i ^k+i^ldxk. Introducing appropriate C°° coordinates, X can 
be written as: 



(3.1) 


with f(x) = 0(||x|p) where x = (xi,.. .,x„). It is said that X has a nilpotent singularity of 
codimension n at 0 if the generic condition d^fldx^(0) 0 is fulfilled. The vector field X 
itself will be often referred to as a nilpotent singularity of codimension n. 


Nilpotent singularities of codimension n are generic in families depending on at least n 
parameters and according with [18, Lemma 2.1] we can state the following result: 


Lemma 3.1. Any n-parametric generic unfolding of a nilpotent singularity of codimension n in R” 
can be written as 


n-1 


k=l 


pi + ykXk + ^? + h{x, p) 


k=2 


_d_ 

dx„' 


(3.2) 


where p = (pi,..., p„) e R”, h{0, p) = 0, (dh/dxi)(0, p) = 0 for i = l,...,n, (d^hj 9x^(0, p) = 0, 
h(x,p) = OiWix, p)f) and h(x, p) = 0(||(v2, ... ,v„)||). 


Remark 3.2. Besides the condition d^fldx^(0) ^ 0 in (3.1), genericity assumptions in Lemma 3.1 
include a transversality condition involving derivatives of the family with respect to parameters. 
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The classical techniques of reduction to normal forms could be used to remove terms in 
the Taylor expansion of h but we do not need to work with simpler expressions. To obtain 
the results provided in the next sections we will have to impose 

d^h 

K=^^(0,0)A0, (3.3) 

0x10x2 

as an additional generic assumption. 


3.2. Rescalings and limit families. Generalizing the techniques used in [14] for dimension 
three, we rescale variables and parameters by means of 

pi = 

pjc = for k = 2,..., n, (3.4) 

Xk = fork = l,...,n. 


with £ > 0 and + ... + = 1, and also multiply the whole family by a factor 1/e. In new 

coordinates and parameters (3.2) can be written as 




(3.5) 


with K as introduced in (3.3) and where y = ( 1 / 1 , ...yn) belongs to an arbitrarily big compact 
inR". 


The first step to understand the dynamics arising in generic unfoldings of n-dimensional 
nilpotent singularities of codimension n is the study of the bifurcation diagram of the limit 
family 


n-1 


k=l 


dyk 


n 

+ (vi + ^ Vkyk + 

k=2 



(3.6) 


obtained by taking £ = 0 in (3.5). Structurally stable behaviours and generic bifurcations in 

(3.6) should persist in (3.5) for £ > 0 small enough. 


If vi > 0 then (3.6) has no equilibrium points. Moreover the function 


L(l/ 1 , ...,yn) = yn- Viyi - V 3 I /2 - ... - Vnyn -1 


is strictly increasing along the orbits and therefore the maximal compact invariant set is 
empty. Hence we only need to pay attention to the case vi < 0. 

On the other hand, up to a change of sign, family (3.6) is invariant under the transformation 
(v,y) (vi,(-!)"■ V 2 ,(-l)"“^V 3 ,...,v„_i,-v„, 


(-1)"1/d(-1)" V2 ,(-1)" ^y3,...,yn-i,-yn), 


(3.7) 


with V = (vi,..., v„). As a first consequence, the study of bifurcations can be reduced to the 
region 

R = {(vi,..., v„) € S''“^ : vi < 0, Vn < Oj. 
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Moreover, since the limit family is invariant under (3.7) up to a change of sign, for 
parameter values on the set 

T = {(vi,...,v„) e : v„_ 2 , = 0 with / = 0,..., L(n - 2)/2J}, 


where [-J denotes the floor function, the correspondent vector fields in the limit family (3.6) 
are time-reversible with respect to the involution 

R : (l/i, 1/2,1/3, • • •, Vn) ^ ((-l)"l/l, (-l)"“^l/2, • • •, Vn-l, -yn)- 

We said that the manifold T' of dimension \ nl2\ -1 is the reversibility set of the n-dimensional 
nilpotent limit family. 

Note that the divergence of the limit family (3.6) takes the constant value v„. Therefore the 
condition v„ = 0 characterizes a subfamily of volume-preserving vector fields. Assuming 
that n is even and defining m = n/2, for parameter values on the reversibility set the limit 
family (3.6) can be written as 

n-l ^ m-1 ^ 

H ^^+1^ + ('^1 + X ^2k+iy2k+i + 3/i)^- (3-8) 

k=l k=l 

In Appendix A we will prove the following result. 


Theorem 3.3. Introducing the new variables q = S ■ {yi, 1 / 3 ,..., yn-iY (ind p = ( 1 / 2 , y^,---,ynf, 
with 

' -V 3 -V 5 ••• -Vn-l 1 ' 

-V5 . • ■ . • ■ 0 

s= : .•• .•• .•• : , 

~1^n-l 

I 1 0 . oj 

the family (3.8) transforms into 

dp dq dq dp' 

where 

1 

Hiq,p) = 2 < Sp,P > +Viq). 

The potential V is defined as 



^ m-1 

2 y, Vlk+lh+l^m 

^ k=l 


m-1 m-1 


LE ^Ik+lbi-m+k+iyi^n 


k=l i=m-k 


2 Lm/2J 
;=i 

\mll\ m-j-1 

L E biqiqm-j L’ll/m, 

;=1 i=j 


where, given bi = 1, 


i-l 

bi = V 2 (m-i+/;)+i be for i = 2 ,..., m. 

e=i 
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Figure 2. The reversibility curve T is split into several arcs attending to the type of eigenvalues of 
the linear part of (4.2) at p— 


4. Nilpotent singularity of codimension 4 IN 

We will prove that in any generic unfolding of a nilpotent singularity of codimension four 
in R"^ there exists a bifurcation hypersurface of homoclinic connections to bifocus equilibria. 

Along this section we will take n = 4 in all the general expressions introduced in §3. It fol¬ 
lows from Lemma 3.1 that any generic unfolding of the nilpotent singularity of codimension 
four in R^ can be written as in (3.2). After applying the rescaling (3.4) we get 

^ ^ ^ ^ ^ ^ ^ ^ + 0(e^))^, (4.1) 

withv = (vi,V 2 ,V 3 ,V 4 ) e and £ > 0. 

As mentioned in §3.2 the first step to understand the dynamics arising in (4.1) is the study 
of the limit family 

obtained from (4.1) taking £ = 0. As argued in §3.2 one only need to pay attention to 
parameters in the region R = {(vi,V 2 ,V 3 ,V 4 ) e : vi < 0, V 4 < 0}. When v € R, vector fields 
in the limit family (4.2) have equilibrium points p+ = (+ yj-vi, 0,0,0) with characteristic 
equations 

r^ - v/pr'^ - v^r^ - V2r + 2 ^/-Vl = 0. (4.3) 

Local bifurcations arising in the family were discussed in [19]. 

For parameters on the reversibility curve 7“ = |(vi,V 2 ,V 3 ,V 4 ) e : V 2 = V 4 = o| withvj < 0, 
the characteristic equations reduces to r^ - + 2 ^/-vi = 0. It follows that the linear part at 

p+ always have a pair of real eigenvalues and a pair of complex eigenvalues with non-zero 
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real part. Local behaviour at p- is richer and it is depicted in Figure 2. Note that we only 
have to pay attention to = 1 with vi < 0. It easily follows that the linear part at p- has 

• a double zero eigenvalue and eigenvalues +1 at BT = (0,0, 1,0), 

• a double zero eigenvalue and a pair of pure imaginary eigenvalues at HDZ = 
( 0 , 0 ,-L 0 ), 

• two double real eigenvalues at BD = (vi, 0 ,V 3 , 0 ) with v^-8 V-vi = 0 and 

V3 > 0, 

• two double pure imaginary eigenvalues ±i{—V3/2)^/^ at HH = (vi, 0 ,V 3 ,0) with - 
8 ^/-Vl = 0 and V 3 < 0 , 

• four non-zero real eigenvalues ±Aj^, with k = 1,2 for parameters along the open arc 
Sli between BD and BT, 

• four complex eigenvalues with non-zero real part p + coi and —p ± coi for parameters 
along the open arc between BD and HH, 

• four pure imaginary eigenvalues with k = 1,2, for parameters along the open 
arc ‘TL'TT between HH and HDZ. 


From the analysis of the linear part at the equilibrium points it follows that a bifocus is 
only possible at p-. In order to show the existence of bifocal homoclinic bifurcations in the 
unfolding of the nilpotent singularity of codimension four in we must study the existence 
of homoclinic orbits to p- for parameter values along 7“. 

To study the family (4.1) close to the reversibility curve 7“ with vi < 0 it is more convenient 
to use a directional version of the rescaling (3.4) taking vi = -1 and (v 2 , V 3 , V 4 ) = (v 2 , V 3 , V 4 ) e 
R^ to get 

1 / 2 ^ + 1 / 3 ^ + y 4 ^ + ( - 1 + V 2 y 2 + V 3 y 3 + V 4 y 4 + y\ + e^yiyi + 0(£^))^. (4.4) 

The equilibrium points when £ = 0 are given by q± = (+1,0,0,0). Note that in fact q± are 
the only equilibrium points even for £ > 0 because in (3.2) h{x,p) = 0 (||(v 2 ,.. ■,x„)\\) and 
this property is preserved by the rescaling. In order to compare with equations already 
considered in the literature we translate q- to the origin applying the change of coordinates 

xi = {yi + l)/2, V2 = 1/2/2^/^, X3 = 1/3/2^^^, V4 = ^4/2^/"^, 


to (4.4) and multiplying by the factor 2^/^ to obtain 


V2^ + V3^ + + ( - Vi + 7]2^2 + 73^3 + 74^4 + xj + eKXiX2 + 0(£^))^ (4.5) 


dX3 


with ?]2 = 2“^/^(v2 - £k), 7]3 = 2“^/^V3, 7]4 = 2“^^'^V4 and £ = 2^/^£. The equilibrium point q- 
in (4.4) corresponds to the equilibrium point of (4.5) at the origin. The limit subfamily for 
7]2 = P 4 = £ = 0 is now given as 

d d d ^ ^ 

X 2 3 - 1 - X 3 3-h X 4 3-h 

dxi dX2 0X3 




(4.6) 
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Writing u = Xi, (4.6) is equivalent to the fourth order differential equation 

+ Pu"{t) + uit) - ii{tf = 0, (4.7) 

with P = —ri 3 . As already mentioned in the introduction, the above equation has been 
extensively studied in the literature. 

In [1] authors prove that (4.7) can be written as a hamiltonian system (as we have stated in 
Theorem 3.3 for a more general case) satisfying the hypothesis required in [24, Theorem 2] 
to conclude that, for each P < -2, there exists an even solution u with u{t) —> 0 when t —> +oo 
satisfying that u > 0, u' < 0 and (P/2 )m' + u"' > 0 on (0, oo). They also prove that for all 
P < -2 any such even solution is unique. From [3] it follows that this unique homoclinic 
orbit is transversal for the restriction to the level surface of the hamiltonian function which 
contains it and, consequently, it is non degenerate in the sense of Definition 2.14. Moreover, 
again in [1], the persistence of such homoclmic solutions is argued for P > —2 but close 
enough to -2. Variational methods used in [8] allow to prove that at least one homoclinic 
solution exists for P < 2. On the other hand, in [3, Section 2] authors check all hypothesis 
required in [12, Theorem 4.4] to conclude that a Belyakov-Devaney bifurcation takes place 
at P = -2. It consists in the emerging from the primary homoclinic solution and for each 
n e IN of a finite number of n-modal secondary homoclinics (or n-pulses) which cut n times 
a section transversal to the primary homoclinic orbit [13, 5, 7]. Heuristic arguments in [3], 
supported by numerical results, show that the non-degenerate n-modal homoclinic orbits 
arising at P = -2 become in degenerate orbits and disappear gradually when P varies from 
P = -2 to P = 2 through a cascade of coalescences and bifurcations. In particular, it is known 
from [27] that for P close to P = 2 there exist at least two even homoclinic solutions and from 
the numerical results it seems that no other homoclinic orbits reaches P = 2. 

All the above results about the existence of homoclinic solutions of (4.7) can be directly 
translated to family (4.6) and also to the reversible subfamily of (4.2) obtained restricting to 
parameter values along the previously defined reversibility curve 7“. For the later case we 
can conclude that (see Figure 2) 

• for parameter values along U [BD] U S'R there exists a symmetric homoclinic orbit 
at p- which is unique and non degenerate along [BD] U S'R, 

• BD is a Belyakov-Devaney bifurcation point, 

• numerical continuation shows that the non degenerate n-modal homoclinic orbits 
arising at BD become in degenerate orbits and disappear gradually when parameters 
move along D'p' in the direction of HH. Close to that point only two symmetric 
homoclinic orbits persist. 

To study the persistence of homoclmic orbits we will consider (4.1) as an unfolding of 
the Belyakov-Devaney bifurcation point BD. As already mentioned it is better to work with 
expression (4.5) for the rescaled unfolding. With respect to parameters {r] 2 , pa, ?] 4 , e) the point 
BD corresponds to (0,2,0,0). Note that (4.5) can be written as 


x' = f{x) +g{A,x), 


( 4 . 8 ) 


HETEROCLINIC CYCLES ARISING IN GENERIC UNEOLDINGS OE NILPOTENT SINGULARITIES 


19 


where A = (Ai, A 2 , A 3 , A 4 ) = ( 7 ] 2 , 7]3 - 2, ?] 4 , e), 


fix) = iX2,X3,X4, -Xi + 2X3 + xj) 


and 


giA,x) = (0,0,0, A 1 X 2 + A 2 X 3 + A 3 Z 4 + A 4 KX 3 X 2 + 0(AI)). 

As already mentioned, q± are the only equilibrium points of (4.4) for all £ > 0 and hence 
^(A, 0) = 0 for all A. Observe that only bifurcations occurring inside the region of parameters 
with A 4 > 0 will be observed in the unfolding of the singularity Family (4.8) fulfills all the 
hypothesis imposed to (2.1). In particular, x' = fix) satisfies the following: 

(BDl): It has a first integral 



(BD2): It is time reversible with respect to 


R : (Zi,Z2,X3,X4) HA (Zi,-X2,X3,-Z4). 


(BD3): The origin is a hyperbolic equilibrium point at which the linear part has a pair of 
double real eigenvalues + 1 . 

(BD4): Accordingto[l],thereexistsanondegeneratehomoclmicorbit 7 = [pit) = (pi(f),p 2 (f)/F 3 (f)/F 4 (f)) : 
f € R} to the origin such that pi(f) and p 3 (f) are even functions and p 2 it) and p^it) are 
odd functions and, moreover, pi > 0 ,p 2 < 0 and ^4 - p 2 > 0 on ( 0 , 00 ). 

(BD5): According to Proposition 2.13, since y is non degenerate, both the variational equation 
z' = Dfipit))z and its adjoint z' = -Dfipit))*z has a unique non trivial linearly 
independent bounded solution. The function (pit) = fipit)) is a bounded solution of 
the variational equation and 


(pit) = VHipit)) = ipiit) - piit)^,p4it) - 2p2it),-p3it),p2it)) 


is a bounded solution of adjoint equation. 

Finally, let us consider the bifurcation equation for homoclinic solutions ^°°iA) = 0, with 
: A —> R and A c R^ a neighbourhood of the origin, as introduced in Lemma 2.16. It 
follows from Theorem 2.17 that under the generic condition 




where 



then (4.8) has homoclinic orbits (continuation of y) for parameters on a hypersurface ‘TTom 
with tangent subspace at A = 0 given by 


■^AjAi + <^A2A2 + '^A3A3 + fAiA^ = 0. 


( 4 . 9 ) 
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Note that 


i 

II 

r^oo 

1 plit) dt, 

/ —CO 

^A2 

<^A 3 

^00 

1 F 2 (f)F 4 (f) dt, 

/ —00 

II 


X oo 

P2(f)P3(0 dt, 

CO 


f 


Kpiit)plit) dt. 


Clearly + 0. Since ^ 2^3 is an odd function = 0- Integrating by parts one gets 
■^A3 = - P 3 (t)^ dt A 0. Finally, since pi is a positive function, we also get that A 0. 

Therefore the tangent subspace (4.9) intersects A 4 = 0 transversely Consequently ‘TToin also 
meets A 4 = 0 transversely 


Now we have to study the eigenvalues at the equilibrium point in order to determine 
which types of homoclinic orbits can be unfolded by the singularity Since for A = 0 the 
linear part at v = 0 has a pair of double real eigenvalues +1 and dimW®( 0 ) = dimW“( 0 ) = 2 , 
for all A small enough, then we can expect three different types of equilibrium: a focus-focus 
(bifocus), a node-node or a focus-node. It easily follows that the characteristic polynomial 
at V = 0 is given by 

Q(r, A) = - D{A)r^ - C{A)j^ - B{A)r - A(A), 


with 

A(A) = -1 + 0(A2) B(A) = Ai+0(A2) 

C(A) = 2 + A 2 + 0{Al) D(A) = As + 0(A2). 

The condition for an improper node is given by the discriminant equations 

Q(r,A) = 0, ^(r,A) = 0. 

Note that (r. A) = (+1,0) are both solutions of the discriminant equations. Now it follows 
from a straightforward application of the Implicit Function Theorem that there exist two hy¬ 
persurfaces Tr and D"*" through the origin in the parameter space such that, for parameter 
values on T)~ (resp. the equilibrium point at the origin has a double negative (resp. pos¬ 
itive) real eigenvalue. Moreover the respective tangent subspaces at A = 0 are Ai - A2 + A3 = 0 
and Ai + A 2 + As = 0. Let N.^om = ^Aj/ = ih-h h 0) and N£)+ = (1,1,1,0) 

be the normal vectors to the tangent spaces of 'Horn, D~ and at A = 0, respectively. 
Moreover denote N^^=o = (0/0/0/ !)• Since rank{N<f{o^,N£)-,N;{^=o) = 3/ there exists a surface 
Hom~ = Horn n !D~ transverse to A4 = 0 of homoclinic orbits to an equilibrium point with a 
double negative real eigenvalue. Moreover, since rank{Nf{o^,N£)+,N;{^=o) = 3, there exists a 
surface Hom'^ = Horn n D* transverse to A 4 = 0 of homoclinic orbits to an equilibrium point 
with a double positive real eigenvalue. On the other hand rank(N.^om/hi£)-,N£)+,NA 4 =o) = 4 
if and only if <^Ai “ <^A 3 ^ 0. But, taking into account that p 2 and p 4 are odd functions and also 
that p 2 <0 and p^ - p 2 > 0 on ( 0 , 00 ) it follows that 

X OO 

P2it)iP2it) - P4it)) dt > 0 . 

CO 

Hence we can conclude that rank(N.^om/ hi©-, N£)+, Na 4 =o) = 4. Therefore there exists a curve 
Horn- = Horn n !D~ n transverse to A 4 = 0 of homoclinic orbits to an equilibrium point 
with a pair of double real eigenvalues one positive and the other negative. 
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Summarizing, we have proved the following result 

Theorem 4.1. In a neighbourhood of A = 0 there exists a bifurcation hypersurface "Horn correspond¬ 
ing to parameter values for which (4.8) has homoclinic orbits to the origin. Moreover there exist 
two bifurcation surfaces 9Iom~ and 9Iom^ contained in "Horn corresponding to parameter values for 
which the origin has a double negative and positive, respectively, real eigenvalue. The surfaces 9Iom^ 
and 'T(om~ intersect transversely along a curve "Horn- corresponding to parameter values for which 
the origin has a pair of double real eigenvalues {ri, r 2 } with ri <0 < r 2 . T-iom~ U 9Iom'^ splits "Horn 
into four regions: 

(i) 'Homff: homoclinic orbits to a focus-focus equilibrium (bifocus case), 

(ii) ‘TYojwjv+f-- homoclinic orbits to a (repelling) node-(attracting) focus equilibrium, 

(Hi) ‘7Yo?wf+n-- homoclinic orbits to a (repelling) focus-(attracting) node equilibrium, 

(iv) ‘TYojwjvn- homoclinic orbits to a node-node equilibrium. 

All bifurcations are transverse to A 4 = 0. 

Since all bifurcations are transverse to A 4 = 0 they are also present in the unfolding of the 
nilpotent singularity of codimension four. Particularly, Theorem B follows as a corollary of 
Theorem 4.1. 

Remark 4.2. Recall that the bifurcation point A = 0 in (4.8) corresponds to the Belyakov-Devaney 
bifurcation point BD in (4.1). In particular, the hypersurface of parameters corresponding to homo¬ 
clinic orbits to a node-node equilibrium point in family (4.1) cuts e = 0 along the curve SR. This 
follows by restricting the bifurcation analysis to A4 = 0. 

5. Nilpotent singularity of codimension 3 in R^. 

We will prove that in any generic unfolding of the nilpotent singularity of codimension 
three in R^ there exists a one-side bifurcation curve of topological Bykov cycles. 

Along this section we will take n = 3 in all the general expressions introduced in §3. It fol¬ 
lows from Lemma 3.1 that any generic unfolding of the nilpotent singularity of codimension 
three in R^ can be written as in (3.2). After applying the rescaling (3.4) we get 

1/2^ + 1/3^ + (vi + V21/2 + V31/3 + y\ + £Kyiy2 + 0 (£ 2 ))^. ( 5 . 1 ) 

withv = (vi,V 2 ,V 3 ) € and £ > 0. 

As mentioned in §3.2 the first step to understand the dynamics arising in (5.1) is the study 
of the limit family 

The above family has already been treated in the literature and discussions about several 
aspects of the dynamics can be seen in [14, 15, 16, 29] and references there included. As 
in the case of the nilpotent singularity of codimension four in R"^ we are interested in the 
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dynamics close to the reversibility curve T = {(vi,V 2 ,V 3 ) e : V 3 = 0}. Particularly, we will 
pay attention to parameters with vi < 0 and V 2 < 0 . 

To study family (5.1) close to the reversibility curve it is more convenient to use a directional 
version of the rescaling (3.4) taking V 2 = -1 and (vi, V 3 ) = (vi, V 3 ) € to get 


y2^ + + (vi -1/2 + vm + y? + £Kyiy2 + o(£^)) ^ 


^ys' 


Moreover, in order to use results already present in the literature, we introduce new variables 
(^ 1 /^ 2 /^ 3 ) = -2(yi, y 2 , y 3 ) and write - 2 vi = when vi < 0 to obtain 


S d ( j _ 7 n ^ 

X2 — + X3-^ + [r - X2 + V3X3 - -X^ - 2 eKXiX 2 + 0 (£^)j- 


(5.2) 


'dx3 dx2 ' 2 ^ '9x3' 

In the above expression, taking the limit case e = 0 and also V 3 = 0 we get the 1-parameter 
family 

^ d / rj \ ^\ d 


xi 


^x^ 




9X3 


As mentioned in the introduction the above family has been extensively studied in the 
literature. It is commonly referred as the Michelson system and it has the following properties: 


(Ml): For all c > 0 the system has only two equilibrium points Q+ = (+ V2c, 0,0) where the 
characteristic equation of the linear part is given by -I- r + V2c. 

(M2): The eigenvalues at Q+ (resp. Q-) are A and -p ± ico (resp. -A and p ± ico) with A > 0, 
p > 0 and cu A 0. Therefore dimW®(Q+) = dimW“(Q-) = 2. Moreover, since the 
Michelson system has zero divergence, A - 2p = 0 and hence p < A. 

(M3): It is time-reversible with respect to the involution 


R : (Vi,V 2 ,V 3 ) (-Xi,V 2 ,-X 3 ). 


(M4): It follows from [31] that the Michelson system has a solution p{t) given by 

pi(f) = a(-9tanh|Sf -l- lltanh^/lf), p 2 {i) = Psit) = p”(f) (5.3) 

with n = 15 -y/ll/19^ and = Vll/19/2, when c = = V2a. It parametrizes an orbit 

Ti along which two branches of the 1-dimensional invariant manifolds coincide. The 
orbit Ti is invariant by the reversibility, that is, 

pi-t) = (PlH),P2(-t),p3(-t)) = Rp{t) = i-piit),p2it),-p3it))- 

Hence pi and p 3 are odd functions and p 2 is an even function. 

(M5): It follows from [29] that when c = Ck there also exists an orbit r 2 = W®(Q+) n W“(Q_). 
Moreover the intersection is topologically transversal. 


Putting together (M4) and (M5) it follows that 


(M 6 ): When c = the Michelson system has a topological Bykov cycle. 
Let us recall the notion of topological Bykov cycle. 
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Definition 5.1. Consider a vector field in with two hyperbolic equilibrium points p± satisfying 
dimW'^(p+) = dimW'“(p_) = 2. Any heteroclinic cycle consisting of two heteroclinic orbits Fi c 
W“(p+) n W^{p-) and 12 c W^{p+) n W“{p-) between the equilibrium points p± such that the 
intersection along T 2 is transversal is called a T-point. When at p+ the linear part has complex 
eigenvalues the T-point is called a Bykov cycle. In both cases, if the intersection along T 2 is only 
topologically transversal we refer to a topological T-point or a topological Bykov cycle. 

Remark 5.2. Let us recall the notion of topological transversality. Let Dbea 2-dimensional open disk 
transversal to the flow at some point T e T 2 . Let S and U be the connected components ofW^{P-) n D 
and W®(P+) n D, respectively, containing the point T. If the diameter ofD is small enough both S 
and U split D into two connected components. The connection T 2 is said topologically transversal 
if the two connected components ofS \ {T} belong to different connected components ofD \ U. 

To study the persistence of the Bykov cycle we will consider (5.2) as an unfolding of the 
Michelson system at the Kuramoto point c = Ck introduced in (M4). Note that (5.2) can be 
written as 

x'= f{x) +g{A,x), (5.4) 

where A = (Ai, A 2 , A 3 ) = (c^ - cl,V 3 ,e), 

7 1 7 

fix) = iX 2 ,X 3 ,C^-X 2 - -X^) 

and 

giA,x) = (0,0, Ai + A 2 X 3 - 2 A 3 Z 1 X 2 + OiAl)). 

Bifurcations occurring inside the region of parameters with A 3 > 0 will be observed in 
the unfolding of the singularity. Family (5.4) fulfills all the hypothesis imposed to (2.1). 
Particularly x' = fix) satisfies all properties from (Ml) to (M 6 ). 

The heteroclinic orbit Fi is non degenerate of codimension two in the sense of Defini¬ 
tion 2.14. Hence the variational equation z'(f) = D/(p(f))z(f) has a unique (up to mul¬ 
tiplicative constants) bounded solution fipit)) whereas the adjoint variational equation 
w' it) = -D fipit))* wit) has a pair of linearly independent bounded solutions, again according 
to Proposition 2.13. Let (pit) = i(piit),(p 2 it),(p 3 it)) and (pit) = i(piit),(p 2 it),(p 3 it)) two of such 
solutions. Since (pit) A (pit) is a bounded solution of the variational equation it follows that 
the plane determined by (pit) and (pit) is orthogonal to fipit)) for all values of t. Therefore, all 
solutions of the adjoint variational equations with initial conditions on /(pjO))-*- are bounded 
solutions. 

On the other hand, it easily follows that w'it) = -Dfipit))*wit) is invariant under the 
involutions iw 3 ,W 2 ,W 3 ) i-> i-W 3 ,W 2 ,-W 3 ) and iwi,W 2 ,W 3 ) i-> iwi,-W 2 ,W 3 ) and the time 
reverse t i-> -t. Therefore, taking (p(0) = (0,-1,0) and (piO) = (1 - c^/p 2 ( 0 ), 0 ,1) we can 
conclude that (pit) and (pit) are bounded solutions of the adjoint variational equation and 
that (pi, (p 3 and (p 2 are odd functions whereas (p 2 , (pi and (p 3 are even functions. 

Consider now the bifurcation equation for homodinic solutions ^°°iA) = 0, with = 
if 1 ,^ 2 ) : A ^ R^and Ac a neighbourhood of the origin, as introduced in Lemma 2.16. 
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Figure 3. Sketch showing the bifurcation curve fiet to topological Bykov cycles in family (5.4) and 
the spirals of bifurcation to homoclinic orbits 'Hom+ and 9dom-for Q+ and Q_, respectively. Note 
that when A 3 = 0, the existence of'Horn- follows from the symmetry with respect to (3.7). 


It follows from Theorem 2.17 that, under the generic condition rankD,\^°°(0) = 2, where 

with 

( 0 ) = 1 

^l,Ai 

^2,Ai 

^ 1 ,A 2 ^ 1 ,A 

^ 2 ,A 2 ^2 ,a 

:) 

<^1,31 = 

1 (p3{t) dt, 

/ —CO 


iy. = J 

f f3{t)dt, 

— 00 

<^1,32 = 

^00 

1 (P3{t)p3(t)dt, 

/ —CO 


= J 

^00 

f f3it)p3{t) dt, (5.5) 

— 00 

^1,A3 ~ 

^00 

1 -2K(p3{t)pi{t)p2{t)dt, 

/ —CO 


i2.A, = J 

^00 

I -2Kf3(f)pi(f)p2(t)df, 

— CO 


then (5.4) has heteroclinic orbits (continuation of Fi) for parameters on a bifurcation curve 
f-let with tangent subspace at A = 0 given by the intersection of the planes 


+ <^1,A2^2 + <^1,A3^3 - 0 
‘^2,AiAi + <^2,A2^2 + <^2,A3^3 = 0. 

From the parities of pi, p 2 , P 3 , (p 3 and 1/^3 it follows that 

■^l.Ai = = <^2,A3 = 0. 

On the other hand, in Appendix B we will show that 

^1,A2 ^ 0/ ^1,A3 ^ 0, (^2,Ai t 0. 

Hence, it is a straightforward application of the Implicit Function Theorem that, indeed, 
there exists a bifurcation curve ffet of heteroclinic connections along the one dimensional 
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invariant manifolds. Moreover it easily follows that the tangent space at A = 0 is generated 
by a vector (0 ,-^i,a 3 /^i,A 2 / 1 ) hence "Het intersects A 3 = 0 transversely Since T 2 is a 
topologically transverse intersection, TYef is a bifurcation curve of topological Bykov cycles. 
This concludes the proof of Theorem C which was stated in the introduction. 

As already mentioned there is a ShiTnikov bifurcation surface 'Hom+ shaped as a scroll 
around TTef (see Figure 3) corresponding to parameter values for which the system has a 
ShiTnikov homoclinic orbit to Q+. Note that the ShiTnikov condition is open and hence 
it follows from (M2). Moreover since the trace of the linear part at Q+ is given by A 2 the 
dissipative condition is also satisfied in 'Hom+ n {A € : A 2 < 0}. Hence strange attractors 

exist for parameter values on positive Lebesgue measure set. 

Remark 5.3. In [16] the existence of subsidiary Bykov cycles in the Michelson system for values 
of c close to Ck is also discussed. Moreover, in [35] the accumulation of Bykov cycles when c —> 0 
is also argued. In all cases the heteroclinic connections have the symmetry properties that we have 
just used. It should be possible to extend our result to conclude that in (5.2), and consequently in 
(5.1), there exist more bifurcation curves to Bykov cycles and particularly an infinite sequence of such 
type of bifurcation curves. Nevertheless, the generic conditions on the bifurcation equation need to be 
checked. 
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Appendix A. Proof of Theorem 3.3 


We will use the following technical result: 

Lemma A.l. Given a symmetric upper anti-triangular matrix 

am U-m-l ■ ■ ■ U.2 1 


A = 


^m-l 


fl2 

1 0 


0 


A ^ is a lower anti-triangular symmetric matrix 


0 


A 


-1 


where, given hi = \, 


for i = 2,... ,m. 


... 0 1 

; ••• ••• b2 

0 • ■ hm-l 

1 &2 ■•■ hm-l hm 

i-1 

hi = - 'Y'^ai-f+ibf 
e=i 


Proof. Let P be an anti-diagonal matrix with all entries equal to 1. Hence 


^ 1 0 ... 

fl2 


L = PA = 


... 0 


Om-l 

^ ^m-1 


'. 0 
fl2 1 


is a lower triangular matrix. Therefore, L~^ = {hi j) is also a lower triangular matrix and hence 
A~^ = L~^P~^ = L~^P is a lower anti-triangular matrix. In fact, using the well know formulas 
for the calculation of the inverse of a triangular matrix, it follows that, for all / = 1 ,.. . ,m 


bij = 0 for all I = 1,..., / - 1, 
i-1 

bi,j = -^ ai-{+ihey for alli = ; -i- 1 ,..., m. 
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On the other hand, bij = &;+i,;+i for all f = ; + 1,..., m - 1. Indeed it is clear for i = j + 1. 
For i = j + 2,... ,m - 1 we can argue by induction. 

Finally, by defining bi = bi^i for all f = 1,..., m, and calculating A~^ = L~^P the proof is 
finished. □ 


It follows from Lemma A.l that 


—1 


1 

bi 

bm-l 


1 &2 • • • bjfi-i bfj 


where, defining bi = 1, 


i-i 


bi = ^ V 2 (m-i+{)+ib{ for f = 2,..., m. 


e=i 


or equivalently 


m-j 


tn-1 


bm-j+i = ^ V 2 (j+t-i)+ibi = ^ V 2 j,+ibk-j+i for 7 = 1,..., - 1. 

t=l k=j 

Writing family (3.8) in the new variables we get 

m-1 m m-1 

\ a / 1 \ o 


(A.1) 


Yj ^i-m+k+l(h) ^ + (^1 + 2], W/c+iPJc + qi) 


dpk 




k=l i=m-k ^ k=l 

To obtain a function V{q) such that p = —W{q) we need -dVjdqi = p, for all i = . ,m. In 

particular 

. 2 

- + 2^ Vlk+lpk + 'Im/ 


dqn 


k=l 


and therefore 


m-1 ^ m-1 

— V{q) = Viqjfi + ^lk+l{^'^k+iqm 2j bi-m+k+iqi^m^ 


k=l i=m-k 

1 3 

+ ^qm + (Pm-liqi, ■ ■ ■ ,qm-l)- 

From the identity -dVIdqm-i = pm-i and taking into account the equation (A.l) we get 

-1 m m-1 m-1 

o(pm-l 


dqm- 


1 


i=l 


k=l 


and therefore 


J m-2 

Pm-liqi/ ■ ■ ■ / qm-l) ~ '^m-iqffi_\ biqiqm-l "t (Pm-liqi/ ■ ■ ■ / qm-i)- 

i=l 
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Since -dVIdqm-i = pm-2, a similar computation leads to 

-v m m-1 

OCpm-l 


m-2 


dqm-2 


^2k+l^k-l^tn ^m-2qm-l — 


i=2 


k=2 


i=2 


and hence 


^ m-3 

(Pm-2iqi/ ■ ■ ■ / qm- 2 ) — '^m-3^m-2 + (Pm-si^lr ■ ■ ■ / 

i=2 


A recursive argument provides 


m m-1 m-1 m-] 

"tT ~ / , ^i-j+l^i ~ / , T^2J:+1^A:-;+1 “ / , ^i-j+l^i — / , ^i-j+l^i/ 


m-1 


i=m-j+l 




for all 7 = 1 ,..., Lwi/ 2 J and consequently. 


(Pm-ji^l/ ■ ■ ■ / ^m-j) ~ '^m-2i+l^jji_j 

m-j-1 

+ biqiqjji-j + cpm-j-iiqi/ ■ ■ ■ ,qm-j-i) 

i=i 

where for j = lm/2j the function (pm-Lm/ 2 j-i is constant. Therefore we get a function V{q) 
with 

m-1 m-1 ^ 

-V(q) =Viqm + '^V2k+l\^-h+iqm+ ^ bi^m+k+l^i^mj + ^qh 


k=l 


i=m-k 


Lm/2J 


m-j-1 


+ 


{^2j^m-2i+l^m-i biqiqm-j^ + (Pm-\ml2\-l> 

j=l i=j 

such that p = —W{q). This concludes the proof of Theorem 3.3. 
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Appendix B. Estimation of ^i,A2/ and (^2,Ai 
Let us write the adjoint variational equation w'(t) = -Df{p{t))*w{t) as 

= -ZViit) + Wsit) 
w'^{t) = -W2(t) 


(B.l) 


where zv = (zvi,W 2 ,'iV 3 ), pi(t) = a(-9tanh(|3t) + 11 tanh^(|3t)), with a = 15 Vll/193 and = 
Vll/19/2, and p = ipi,p 2 ,P 3 ), withp2 = p[ and ps = p^'• Writing s{t) = tanh(|3t) it follows that 

piit) = njSjl - s{t)^){-9 + 33s(f)^) and p3(f) = - s(f)^)(84s(f) - 132s(f)^). 

Our goal is to show that ^2,Ai' given in (5.5), are different from zero. Note 

that from the parities of pi, p2, ps, 1/^3 and cp^, it follows that 


^2, 
<^1,A2 


^00 

,,Ai = 2 I 1/^3 (f) dt, 

Jo 

^00 

2 I (p3it)p3it) dt, 

-4k 

Jo 


(B.2) 


(p3it)Pl(t)p2(t) dt. 


We have already argued in §5 that each bounded solution w{t) satisfies an orthogonality 
condition with respect to /(p(f)), that is, 

P2(t)m(t) + P3(t)W2(t) + {cl - p2it) - ipiit)f/2)W3it) = 0. 

with Ck = V2a. Now, introducing vjt) = zv 3 {t) and V 2 {t) = Dj(f), we get equivalently 


A{t)v'{t) = B{t)v{t), 


(B.3) 


withD(f) = (Di(f),D2(f)) and 


A{t) 


1 0 
0 P2{t) 


Bit) 


0 


{pi{t))^l2-cl 


1 

P3{t) 


Since we have chosen bounded solutions (p(t) and ^{t) of (B.l) satisfying the initial conditions 
(p{0) = (0,-1,0) and i/i(0) = (1 - c^/p2(0),0,1), they correspond, with respect to the new 
variables Vi and V 2 , to solutions (p{t) = {(pi{t),(pi{t)) and ipit) = (i/'i(f), i/'2(f)) of (B.3) with 
initial conditions (p(0) = (0,1) and Tp{0) = (1,0), respectively. Therefore, taking into account 
that K t 0, we only have to show that 


^00 r^o 

I i/)i(f)df^0, 

Jo Jo 


(pi(t)p 3 (t)dt ^ 0 and 


f 


(piiOpiiOpiit) dt t 0. 


Note that A{t) is singular at t+ = tanh“^( V3/11)/|S ~ +1.5529. Hence (B.3) must be treated 
as a differential algebraic equation rather than as an ordinary differential equation on the 
interval [0,00). We will provide approximate values of the integrals in (B.2) using numerical 
methods on an interval [0, fo]/ with to > t+ to be fixed later, and providing upper bounds for 
the absolute value of the integrals on [fo, 00). 
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In order to get appropriate upper bounds on [to, oo) we write the equation (B.3) as v'{t) = 
Qv{t) + R{t)v{t) with R{t) = - Q and 

Q = ,i“S('^<‘r'«(‘) = (-3o“/i9 -v^)' 

The eigenvalues of Q are A+ = (- ^/209 + i V2071)38 and 

^ _ I (-V^-zV2(m)/60 (-V^ + zV2(m)/60 J 

is such that Q = P}P~^ where / is the complex canonical form of Q with entries A+ and A- 
along the diagonal. From the Lagrange Formula 

v{t) = Pe^^^~^°^P~^vo + Pe^* r e~^^P~^R(s)v{s) ds, 

Jto 

for the solution v{t) of (B.3) with v{to) = vq, we obtain 

e-^^Mt)\\ < \\P\\e-^^°\\P-^\\\\vo\\ + IIPII f c-“1|P-i||||R(s)||||z;(s)||ds, 

Jto 

where a denotes the real part of A+. For any c > 0 we can choose to such such that ||R(t)|| < e 
for all t > to and hence, applying the Gronwall Lemma we get 

Mt)\\ < ||P||||P“^||||r;olle<“+"^""^“'""><*“^“^ 


Moreover we can assume that to is large enough to have ||R(t)|| strictly decreasing on [to, oo) 
and ||P||||P“^||||R(to)|| |fl|. Therefore we can take e = R(to) and achieve the following upper 

bounds 




- X” 11^(011 

< l|P|ll|P“^lll|i/'(fo)ll XX dt 

< L||i/)(to)||, 


with 


and, analogously. 


-||P||||P~^II 
fl + iiPiiiip-iiiiiR(to)ir 


roo 

(pi(t)p3(i)dt 

poo 

(p(t)pi(t)p2(t)dt 


< L|p3(fo)ll|(p(fo)ll, 

< L|pi(to)p 2 (fo)ll|(p(fo)ll, 


assuming that to is chosen such that |p3(t)| and |pi(t)p2(f)l are decreasing on [to, oo). 


To conclude we have to provide estimations of the integrals on [0, to]. One can check 
that, as required, the value to = 20 is such that ||R(t)||, |p3(t)| and |pi(t)p2(f)l are decreasing 
on [to, oo). Moreover the numerical results show that the values ||i/’(20)|| and ||(p(20)|| provide 
small enough upper bounds for the integrals on [to, oo). As already mentioned (B.3) must 
be treated as a differential algebraic equation. The computing environment MATLAB [36] 
provides codes to deal with such kind of equations. For instance, according to [51], the 
odelSs code is an appropriate one to solve DAEs. We have used such method working 
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TOE 

Estimation of 

ipi(t) dt 

Estimation of 

f<Plit)p3(t)dt 

Estimation of 

r20 

fo (pi(t)pi(t)p2(t) dt 

10-3 

-2.69317886 

3.38993680 

2.14437726 

10-4 

-2.66952044 

3.42255576 

2.18337593 

10-3 

-2.65457297 

3.42416164 

2.19168295 

10-6 

-2.65558754 

3.48165747 

2.27650773 

10 -^ 

-2.65577082 

3.42402265 

2.19092815 

10-3 

-2.65586498 

3.42412150 

2.19145569 

lO-'^ 

-2.65592518 

3.42421346 

2.19175351 

10-10 

-2.65594764 

3.42423492 

2.19185340 

10-11 

-2.65596101 

3.42424493 

2.19188683 

10-12 

-2.65596369 

3.42424803 

2.19190343 

10-43 

-2.65596540 

3.42424892 

2.19190641 


Table I. Estimations of the integrals on [0,20] using the code odelSs provided in MATLAB to deal 
with differential algebraic equations. Note that MATLAB uses two different tolerances, absolute and 
relative, to compute optimal steps. V\fe take both equal, with the values indicated in the first column. 


TOE 

Estimation of 
k||i/^(20)|| 

Estimation of 
IP3(20 )|L||(p(20)|| 

Estimation of 
|pi(20)p2(20)|L||(p(20)|| 

10-13 

4.219110x10-2 

2.103834x10-2 

3.321829x10-2 


Table II. Upper bounds of the integrals on the interval [20, oo]. 


with different tolerances to get the results shown in Table I. Once the numerical solution is 
obtained, the integrals were approximated using the trapezoidal rule with the nodes which 
were generated by the numerical algorithm. With the highest tolerance value, the maximum 
step sizes are 0.01355215 and 0.01323550 for the solutions ip and (p, respectively. On the other 
hand one can get accurate estimations of the upper bounds for the integrals on [to, oo]. The 
results taking the approximates values of the solutions at t = to obtained with the highest 
tolerance are shown in Table II. 

Remark B.l. All the estimations of upper bounds are obtained using the maximum norm. 
It can be easily checked that ||P|| = 2, ||P“^|| = V30/109 + 30/ V2071 and ||R(t)|| = (|pi(t)^/2 - 
cl + 30/191 + \p3it) + VTT7T9|)/|p2(f)|. 

Remark B.2. The code ode23t included in MATLAB also solves differential algebraic equa¬ 
tions. Even ode45 can be used to solve our particular equation. The results does not change 
significatively, namely, the differences are below lO""^. We have also used our own algo¬ 
rithms (a Taylor method of order 25 for the variational equation with small fixed step and 
projection on the stable manifold). Again the differences are below lO”"^. 
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Appendix C. Proof of Lemma 2.16 


As it stated in Proposition 2.7 the variational equation z' = Df{p{t))z has exponential 
dichotomy in [to/°o) and (-oo, tg]- Let 

= X{to)P+X-\to) and I - = I - Xito)P-X-\to) 

be the corresponding projection matrix on the stable space = Tp(f(j)W®(p+) and instable 
space £“ = Tp(tg)W“(p_), respectively. 

Before to give the proof of Lemma 2.16 we need the following preliminar result. 

Lemma C.l. Let b e C^{[to, oo), IR”). Then, z'^(t) is a positively bounded solution of 


z'=Dfip{t))z +bit) 


if and only if 


z+(f) = Xit)X-\tQ)LX>+itQ)z^itQ) 


f 


+ I X(t)X-^(s)j^+(s)b(s)ds 


J r»Oi 

t 


Xit)X-\s){I - ‘^+is))b{s)ds. 


(C.l) 


Proof. Since X(f) is the fundamental matrix of the lineal homogeneous equation z' = D/(p(f))z, 
the solution of the complete linear equation z' = D/(p(f))z + bit) are 


z(t) = X(f)X-i(fo)z(to) + X(t) 


/’x-’ 


is)bis) ds. 


By means of the projection P+ in the exponential dichotomy of the homogeneous equation 
hr [fo/ this solution can be written as 

zit) = X(f)P+X-i(fo)z(fo) + Xm - P+)X-\to)zito) 

(C.2) 


f 


+ X(f) I P+X-\s)bis)ds + Xit) 




P+)X-fs)bis)ds. 


On the other hand, according to the exponential dichotomy ||X(f)P+X ^(s)|| < Ke for 
t > s > to, it follows that 

|X(f)P+X-i(fo)z(fo)l < for t > to, 

|X(f) r P+X-\s)bis) ds\ < r ds for f > fg, 

Jfo Jto 

and thus, the first and third term of (C.2) are bounded for t > to. 

If we assume that z(f) is bounded, necessarily then the sum 


X(f)(l-P+)X-i(fo)z(fg) + X(f) I iI-P+)X-fs)bis)ds 


/' 


r-1/ 


= X(f)[(l - P+)X-i(to)z(fo) + I (f - P+)X-^is)bis) ds] 


/'< 


(C.3) 
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is also bounded. However, the exponential dichotomy ||X(t)(7 - P+)X ^(s)|| < Le for 
s >t >to, implies that 

|X(fo)a - P+)X-i(to)z(fo)l = |X(fo)(f - P+)X-i(f)X(f)(f - P+)X-i(to)z(fo)l 

< - P+)X-\to)z{to)\. 


Hence, 

|X(t)a-P+)X-i(fo)z(fo)l > |X(fo)a-P+)X-i(fo)z(fo)|P-^e^*''-*“> 

is not bounded and since |X(f)(J - P+)X“^(fo)z(fo)l ^ l|X(f)|||(f - P+)X“^(fo)z(fo)l/ it follows that 
the matrix X(f) is not bounded. Therefore, from (C.3) we get that the solution z{t) only can 
be bounded if it holds 


(J-P+)X-i(fo)z(fo) + 


^OO 

dfn 


P+)X-\s)b{s)ds = 0. 


Consequently, from (C.2) we obtain that if z'^{t) is a bounded solution of z' = Df{p{t))z + b{t) 
then 

z+(f) = X(f)X-i(fo)^+(to)z+(fo) 


+ 


f 

Ptg 


X(t)X-^(s)^+(s)b(s)ds 


-r 


X(t)X-^(to)(I - ^+(s))b(s)ds. 


Conversely, to verify that z'''(t) given in (C.l) is a bounded solution of z' = Df(p(t))z + b(t) it 
suffices to see that 


L+b(t) 


Jto 


X(t)X-\s)^+(s)b(s)ds- 


Jf 


X(t)X-^(to)(I - ^+(s))b(s)ds 


is a particular bounded solution of the above complete linear equation. Indeed, 

7 ^OO 

-L^b{t) = Df{p{t))X{t)[ P^X-\s)bis)ds- I {I-P+)X-Hs)b{s)ds] 


+ . 


• X{t)[p+X-\t)b{t) + (J - P+)X-Ht)b{t)\ = Df(pit))L+b(t) + b(t), 

J r>f ^OO 

Ke-‘^(^-%{s)\ds + ds<K + L, 

tn Jt 


for all t > to, where the constant K and L not depend on t. 


□ 


In the same way, a similar result to the negative bounded solutions of the complete linear 
equation is followed. 

Lemma C.2. Let b e C°((-oo, fo],K”). Then, z~(t) is a negative bounded solution of 

z'=Dfipit))z + bit) 

if and only if 

z-(t) = X{t)X-Hto){I - t^-ito))z-ito) 

pf rk (C.4) 

+ X{t)X-\s)‘^-{s)b{s)ds- Xit)X-\s){I-^-is))bis)ds. 

%J—oo %J t 
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Remark C.3. Notice that the functions given in (C.l) and (C.4) are both solutions of the equation 
z' - Df{p{t))z + b{t)for any continuous function b(t) on [^ 0 /°°) (ind (- 00 , ^q], respectively. On the 
other hand, to prove that both solutions are bounded solutions of this equation we need to use that 
b(t) is also bounded. 


Recall that = Tp(tuW®(p+)and£“^ = Tp{t^)W^{p-). Notice that since the (homo)heteroclmic 
orbit Y is non degenerate, Efg = n £“ is one-dimensional. In fact, this space is generated 
by the vector p'(to) = fipito))- Moreover, according to Remark 2.15 the dimension of 


E®* n E“* 


K 




is d = S- - s+ -E 1 where s+ are the stable indexes of p+ and p- respectively. 

Let us define W* the ortogonal complement to Etg in the tangent space of the stable and 
unstable manifolds, respectively, that is E®^ = Ef^J.W+ and E“ = Ef^J.W“. Then 

= span{/(p(fo))} © W+ © W" © (C.5) 


Finally, we take the transversal section to y at p{to) given by 


= p{to) + {fipito))}^ = p{to) + [W+ © W- © e;j 
and a base {wi: i = 1... d} de E*^. Notice that 

Wi{s) = X~^{s)*X{to)*Wi for i = 1,... ,d 


are bounded linearly independent solutions of the adjoint variational equation. 


Proof of Lemma 2.16. The solutions P/v(f) of (2.1) can be written as pA(t) = pit) + z^it) where 
z^it) is a solution of (2.2). Let Yfg = W+ © W“ © E*^. Assuming pA(to) e Etg then ZA(to) e Y^g. 

In order to get that p,\(f) parametrizes an orbit in W®(p+(A)) (resp. W“(p_(A))), the function 
ZA(t) have to be a positively (resp. negatively) bounded solution of (2.2). Now if we assume 
that z-(t) are a pair of positively and negatively bounded solutions of (2.2) respectively, then 
b(-,z-(-),A) € C°(/±,R”) where /+ = [fo/°o) and /_ = (-oo,fo]. Thus, according to Lemma C.l 
and Lemma C.2 it must be met that 

z+(f) = X(f)X-i(fo)^+(io)z+(fo) 

r” 1 (c-6) 

+ Xit)X-\s)‘A^+is)bis,z-is),A)ds- X(f)X-i(s)(l - ,!^+(s))&(s,z±(s),A)ds 

Jto Jt 

and 

z-(f) = X(f)X-i(fo)(J-^-(fo))z-(fo) 

N rk (C.7) 

+ Xit)X-\s)^-is)bis,z-is),A)ds- X(f)X-i(s)(J - ,?^-(s))&(s,z±(s),A)ds. 

%J —00 %Jt 


Conversely, according to Remark C.3, the solutions z'''(f) and z“(f) of the integral equa¬ 
tions (C.6) and (C.7) are both solutions of (2.2), but not necessarily bounded. The existence 
of positively and negatively bounded solutions of (C.6) and (C.7) will be proved as an 
application of Implicit Function Theorem. 
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Ler I]"*" = e and i] = {I - ^-{to))z (to) e Equations (C.6) and (C.7) 


can be written in the form 


z- = ^-{z-,^-,K) 


(C.8) 


where : C°(/+,R”) x W* x —> C°(/+,R”). In order to apply the Implicit Function 
Theorem to the equation z - ^-(z-,;]-. A) = 0 notice first that =7^-(0,0,0) = 0. On the other 
hand, Dz^-(z-,;]-, A) : C|J(/+,R”) —> C|J(/+,R") is the null function for z- = rj- = A = 0. 
Indeed, for any h e C°(/+, R”) it holds that 


DzJi^^{z,i],A)h{t) = r 0(f, 

Jto 

-r 


s)^+(s)Dzb(s,z(s),A)h(s) ds 


0(f,s)(J - £P+{s))Dzb{s,z{s), A)h{s) ds. 

Jt 

Since Dzb(s,0,0) = 0 then (0,0,0) = 0. Similarly it follows that DzJ^~(0,0,0) = 0. 

Therefore, there exists 6+ > 0 such that for every rf- e W* and A e R^ with |7]-|, |A| < 6+ there 
is a unique z-{i]-,A) e C°(/±,R”) so that z-(7]-. A) = J^-(z-(i]-,A),ri-,A) and z-(0,0) = 0. 

Now, consider the condition z“(A, rj~)(to) - z'''(A, rj^)(to) e E*^. Since Ytg = © W“ © E*^, 

we can write 

Z“(']“/ A)(fo) = ?]- + w'^ir]-. A) + g-iif. A) 

where w'^(r]-,A) € and Q-(r]-,A) e E*^. From (C.6) and (C.7), and having into account 
that X(t)X~^(s)^±(s) = ^+(f)X(f)X“^(s), it follows that 


z (?] ,A)(fo) = i] +^-(fo) r X(to)X ^(s)&(s,z (?] ,A)(s),A)ds, 

kJ — CO 


f 


(C.9) 


z^(ri^,^)(to) = ri^ - (I- ^+{to)) I X(to)X \s)b(s,z^(ri^,A)(s),A)ds 


Thus 


w'^irj ,A) + q (?] ,A) = ^-(to) r X(fo)X ^(s)&(s,z (?] ,A)(s),A)ds, 

k) —OO 

J rkOO 

X(fo)X“^(s)&(s,z+(7]+, A)(s), A) ds. 

tn 


(C.IO) 


Recall that z-(0,0) = 0 and hence w-(0, 0) = p~(0, 0) = 0. On the other hand, applying that 
Dzb(t, 0,0) = 0 in (C.IO) we get that 

D^.g^(0,0) = D;^g^(0,0) = 0, 

Drj+w-(0,0) = Dazv-(0,0) = 0 . 


(C.ll) 


The condition z“(A, r]~)(to)-z'^(A, r]'^)ito) c E*^ is equivalent to the system of two equations 
rj- - zv-(ri'^. A) = 0 which we write in the form F(t], A) = 0 where i] = (r]'^, ri~) e x W“. We 
have that F(0,0) = 0 and according to (C.ll) it follows that D,jF(0,0) = L Hence, applying 
again the Implicit Function Theorem we get ?] as a function of A. We conclude that there 
exists 6 > 0 (6 < 6+) such that for any A € R^ with |A| < 6 there is a unique ?]“(A) € W* so 
that ?]“(A) - a;-(7]'^(A), A) = 0 and rj-{0) = 0. 
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P(to) + W- 


W" 



PWo) 


P^ito) 




Figure 4. The figure shows the transversal section = p{to) + ®E*^] at the point 

p(io) to a non degenerate (homo)heteroclinic orbit y = {p{t) : t e ]R}. The curves W® and W“ are, 
respectively, the intersections W®(p+(A)) n and W‘‘{p-{A)) n Ef„. A priori the curve W® does not 
met W". However, the projection along the direction of both curves on p(to) + {]Nf © ]Nf ] have a 
unique transversal intersection point q^. Now, q\ + EJ^ mets, respectively, W® and W“ at pf{to) and 
pj(to)- This two points define the vector The persistence of the connection holds if£,°°{A) = 0 

which provides a set of d = dimEJ^ conditions. 


Now consider the functions z*(t) = z-{r]-{A), A){t). The uniqueness, boundedness and 
regularity respect to A of z*(f), and thus of p*(t) = p(t) +z*(t), are following from the Implicit 
Function Theorem. Also we have that zj = 0 and since for |A| small enough z* is close 
to Zq = 0 we get that sup^gj^ ~ F(0I is arbitrarily small. This means that the orbits 
parameterize by p*(f) are close to the orbit y which is parameterized by p{t). So, together 
the hyperbolicity of the equilibria p±(A) we get that limf^+ooP*(f) = p+(A). Thus, for each 
A with |A| < 6, the solutions p*(f) parameterizes, respectively, orbits in the stable manifold 
W^(p+(A)) and in the unstable manifold W"(p_(A)) such that p*(fo) ^ pito) + ^fo = ^fo 
f,°°{A) = pf(to) - F^(fo) = zf(to) - e E*^. This proves the first item of the lemma. 

In order to prove the second item notice that if ^“(A) = 0 then a (homo)heteroclinic orbit 
of (2.1) is given by 



paraf<fo/ 

p+(t) para t > to- 


where p*(0 are the solutions in the first item. On the other hand, if pxit) is a solution 
parametrising a (homo)heteroclinic connection such that p,\ (to) e with|A| and |p/v(to)-p(to)| 
small enough, its restriction to the intervals /_ = (-oo, tg] and /+ = [to, oo) define a pair of 
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solutions in the assumption of the first item. That is, p*(fo) ^ ^to p^(fo) “ P^(fo) = 
0 e E*^. This makes obvious the reciprocal implication. 

To conclude the proof of the second item notice that 

^“(A) = p-(fo) - ptik) = z-^{to) - zlik) € El 
Since {zvj: i = 1... d] is a base of E*^ = n E^* = + E^-^- then 


nA) = '^<wi,nA)>wi. 

i=l 

Form (C.9) and having into account that < Wi, p- >= 0, it follows that 

^“(A) =< Wi,nA) >=< zvi,^.{to) r X(fo)X-i(s)&(s,z-(s),A)ds > 

%J — oo 


'f 


+ <Wi,(I- ,^+(to)) I X(to)X ^(s)b(s,zUs),A)ds > 




< [^-(to)X(to)X ^(s)rwi,b(s,z.(s),A)ds > 


< [(I - ^+(to))X(to)X ^{s)]*Wi,b{s,zlis), A) ds > . 
Thus, since ^-{k)* : R” —> E'^* and I - ^+{k)* '■ R” —> we get that 

[^-ik)Xik)x-\s)Ywi = x-\syx{kT^-ikrm = wys) 

[{I - ^4k))Xik)x-Hs)Ywi = x-Hsyxikni - ^4kr)m = m{s). 
Substituting in (C.12) we obtain that 


^r(^) 


I" 


f 

»7fn 


(C.12) 


< Wj(s),b(s,z,(s),A) > ds + I < Wi(s),b(s,zUs),A) > ds = 0. 


This concludes the second item and proves Lemma 2.16. 


□ 
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